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REACTING VISCOUS-SHOCK-LAYER SOLUTIONS WITH MULTICOMPONENT 
DIFFUSION AND MASS INJECTION* 

By James N. Moss 
Langley Research Center 

SUMMARY 

This study presents numerical solutions of the viscous-shock- layer equations where 
the chemistry is treated as being either frozen, equilibrium, or nonequilibrium. Also the 
effects of the diffusion model, surface catalysis, and mass injection on surface transport 
and flow parameters are considered. The flow is treated as a mixture of five inert and 
thermally perfect species. The viscous-shock-layer equations are solved by using an 
implicit-difference scheme. 

All calculations are for hyperboloids with included angles of 20° and 45°. The flight 
conditions are those for various altitudes and velocities in the Earth’s atmosphere. Data 
are presented to show the effects of the chemical models; diffusion models; surface cataly- 
sis; and mass injection of air on heat transfer; skin friction; shock standoff distance; 
wall pressure distribution; and tangential velocity, temperature, and species profiles. 

The results show that an equilibrium analysis can substantially overpredict the 
heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. More- 
over, at such conditions surface catalysis significantly influences heat-transfer and 
flow-field properties. If a binary rather than a multicomponent diffusion model is 
assumed, negligible errors in most flow properties result. Quantitative results are pre- 
sented that show the effect of mass injection on flow properties within and downstream of 
the injection region. 


INTRODUCTION 

Analysis of the flow about a hypersonic vehicle must account for the interactions 
that occur between the reacting outer flow and the vehicle surface. This problem can be 
complicated in many cases because of either surface catalytic effects or mass injection. 
Furthermore, assumptions concerning the flow field, chemistry, and diffusion models may 
significantly affect the accuracy of an analysis for such a problem. 

Numerical solutions to the aforementioned problem have been either inviscid- 
boundary-layer solutions (refs. 1 to 5, for example) or viscous-shock- layer solutions 

♦This paper is based in part upon a thesis entitled "Solutions for Reacting and Non- 
reacting Viscous Shock Layers With Multicomponent Diffusion and Mass Injection "sub- 
mitted in partial fulfillment of the requirements for the degree of Doctor of Aerospace 
Engineering, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, 
1972. 



(refs. 6 to 18, for example). The viscous-shock- layer equations as proposed by Davis and 
Fliigge-Lotz (ref. 19) result in one set of equations uniformly valid throughout the shock 
layer. Consequently, the viscous- inviscid interactions are accounted for in a straight- 
forward manner. 

The most recent numerical solutions of the viscous-shock- layer equations are those 
given in references 6 to 10. These solutions are for stagnation and downstream flow and, 
with the exception of reference 9, are not restricted to the thin- shock-layer approximations. 
The viscous-shock-layer solutions presented in references 11 to 18 are for flow along 
the stagnation streamline. For the stagnation and downstream solutions, reference 7 con- 
siders atomic and molecular oxygen species by use of finite-rate chemistry; reference 8 
considers the effects of injecting argon, air, and helium into air where air is treated 
as one species and the chemistry is either frozen or equilibrium; reference 9 considers 
nonequilibrium air as a reacting mixture of seven chemical species with constant but 
arbitrary Prandtl and Lewis numbers; and reference 10 extends the viscous-shock- layer 
solutions for both equilibrium and nonequilibrium chemistry to include all the following 
features: downstream solution capability, mass injection, and multicomponent diffusion. 
Also, the effect of the injection of water and ablation species into equilibrium air are 
considered in reference 10. 

This study presents numerical solutions of the viscous-shock-layer equations for 
air. The chemistry is treated as being either frozen, equilibrium, or nonequilibrium. 
Particular attention is devoted to the effects of chemistry model, surface catalysis, diffu- 
sion model, and mass injection. The flow is treated as a mixture of five thermally perfect 
species (O, Og, N, Ng, and NO). All calculations are for hyperboloids with included angles 
of 20° and 45 . The flight conditions are those for various altitudes and velocities in the 
Earth's atmosphere. Data are presented showing the effects of the chemical models; dif- 
fusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; 
shock stand-off distance; wall pressure distribution; and tangential velocity, temperature, 
and species profiles. The results reported in reference 10 and herein are the first stag- 
nation and downstream viscous- shock- layer solutions where the diffusion is treated as 

SYMBOLS 
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coefficients of polynomial curve fits for thermodynamic properties 
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mass diffusion parameter defined by equation (18) 
mass diffusion parameter defined by equation (22) 
skin-friction coefficient (eq, (45)) 
mass fraction of species i, p^/p 
mass fraction of element £ (eq. (10)) 
frozen specific heat of mixture, C. Cp 

specific heat of species i, Cp y/fcp 
multicomponent diffusion coefficients 

binary diffusion coefficients 

free energy of mixture 

free energy of species i 

total enthalpy of mixture (eq. (9)), H*/U*^ 

enthalpy of mixture, J2 Cqh,. 

i 
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enthalpy of species i, h?/U* 
integers 

diffusion mass flux of species i, J? a.*/ p 
diffusion mass flux of element e (eq, (12) ) 
thermal conductivity of mixture, K*/p C* ^ 
thermal conductivity of species i 
matrix coefficients defined by equation (B16) 
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N Le,i 


N Le,ik 


matrix coefficients (see eq. (B15)) 

backward rate constant 

forward rate constant 

arbitrary constant defined by equation (22) 

free- stream Mach number 

molecular weight of species i 

molecular weight of mixture 

mass- injection rate, m*/p*U* 

sum of reacting species and catalytic bodies 

number of chemical reactions 

number of reacting species 

Lewis number defined by equation (19) 

p*C *D*. 

multicomponent Lewis number, 2 — II 

K* 

Prandtl number defined by equation (A2) 
shock Reynolds number, P*U*a*/M * 
free -stream Reynolds number, P*U*a.*/n* 
Stanton number defined by equation (44) 
coordinate measured normal to the body, n*/a* 
parameter defined by equation (C12e) 
parameter defined by equation (C13e) 
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parameter defined by equation (C12f) 
parameter defined by equation (C13f) 
pressure, p*/0°*Utf) 
partial pressure of species i 
wall heat- transfer rate, q*/0°*U* 3 ) 
universal gas constant 

radius measured from axis of symmetry to a point on the body 
surface, r*/a* 

coordinate measured along the body surface, s*/a* 
temperature, T*^U* 3 

reduced temperature defined by equation (B12) 
time 

free- stream velocity 

velocity component tangent to body surface, u*/U* 
velocity component normal to body surface, v*/U* 
mass rate of formation of species i, w?a*y4* U* 
concentration of species i, moles/volume 
mole fraction of species i 
shock angle defined in figure 1 
stoichiometric coefficients for reactants 
quantity defined by equation (55) 
coefficients in equation (C8) 
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stoichiometric coefficients for products 

quantity defined by equation (56) 

coefficients in equation (C27) 

coefficients in equation (C32) 

mole mass ratio of species j (eq. (49)) 

free- stream ratio of specific heats 

number of atoms of 2 th element in species i 

Reynolds number parameter (eq. (A3)) 

maximum energy of attraction of colliding i molecules 

maximum energy of attraction of colliding unlike molecules 

binary Lewis number 

variable binary Lewis number 

constant binary Lewis number 

transformed n coordinate, n/n g 

body angle defined in figure 1 

body curvature, a* a.* 

viscosity of mixture, 

reference viscosity (eqs. (Al)) 

coordinate measured along the body surface, I = s 

density of mixture, p */p* 



collision diameter of species i 


a*. 

1J 
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ij 

Subscripts: 

i,j,k 

« 

o 

r 

s 
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Superscripts: 



o 

j 


collision diameter for unlike species (eq. (B13)) 
quantity defined by equation (B9) 


quantity defined by equation (Clle) 


collision integral for diffusion 


collision integral for viscosity 


ith, jth, and kth species 

eth elements 

stagnation -point value 

denotes rth reaction (see eqs. (52)) 

shock value 

wall value 

free-stream condition 

gas- solid interface value resulting from mass injection (see fig. 2) 

a constant in equation that defines forward rate constant (eq. (53)) 
a constant in equation that defines backward rate constant (eq. (54)) 
standard state (pure substance at 1 atmosphere of pressure) 
zero for plane flow and one for axisymmetric flow 
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quantity divided by its corresponding shock value 


* dimensional quantity 

* total differential 

" shock oriented velocity components (see fig. 1) 

ANALYSIS 

A detailed development of the viscous- shock- layer analysis is presented because of 
the unique features of this analysis for studying problems such as those considered in this 
study and, also, the potential application of this analysis for studying much higher energy 
flow problems. The conservation equations and their associated boundary conditions are 
presented for the laminar viscous shock layer about an axisymmetric or two-dimensional 
body. These equations and boundary conditions account for mass injection and treat the 
diffusion as either multicomponent or binary. Equations convenient for solving frozen, 
equilibrium, and nonequilibrium flow are presented. Nondimensionalizing quantities are 
given, along with the transformation used to facilitate the numerical solution of the 
equations. 


Shock- Layer Equations 

The conservation equations that describe a reacting multicomponent gas mixture 
can be found in the literature. (For example, see ref, 20 or 21.) The viscous- shock- 
layer equations are obtained from the conservation equations by using the same procedure 
as given in reference 6. First, the conservation equations are written in the body- 
oriented coordinate system shown in figure 1. Then these equations are nondimensional- 
ized in each of two flow regions with variables which are of order one. That is, the equa- 
tions are nondimensionalized in the region near the body surface (boundary layer for large 
Reynolds numbers), and then the same set of equations is nondimensionalized in the essen- 
tially inviscid region outside the boundary layer. Terms in each of the two resulting sets 
of equations are retained up to second order in the Reynolds number parameter e . By 
combining these two sets of equations so that terms up to second order in both the inner 
and outer regions are retained, a set of equations uniformly valid to second order in the 
entire shock layer is obtained. 

The terms used to nondimensionalize the viscous -shock-layer equations are given in 
appendix A. The nondimensional viscous- shock- layer equations for a chemically reacting 
multicomponent gas mixture can be written as follows: 
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Species continuity: 


( u aC i * C i\ • 

p _ i+v — = w. - 

\ 1 + nK as an/ 


(1 + nK) (r + n cos 9 )•* 


[(1 + nK) (r + n cos e^jjj (5) 


9 



and state: 


P = 


PT| 


R* 


M*C* 

v P>°° 


( 6 ) 


This set of equations has a hyperbolic-parabolic nature (ref. 6), where the hyper- 
bolic nature comes from the normal momentum equation. If the thin- shock- layer approxi- 
mation is made, the normal momentum equation becomes 

ap _ pu 2 x 
an 1 + me 


When equation (3) is replaced with equation (7), the resulting set of equations is parabolic. 
Consequently, the equations can be solved by using numerical methods similar to those 
used in solving boundary- layer problems. After an initial iteration using equation (7), the 
final flow-field solution is obtained by replacing equation (7) with equation (3); thus the 
thin- shock- layer approximation is removed. 

Equations (1) to (7) are the governing viscous-shock-layer equations. This set of 
equations is convenient to apply when the flow-field chemistry is assumed to be either 
nonequilibrium or frozen. However, when the flow field is assumed to be in chemical 
equilibrium, the production terms w^ that appear in the energy and species continuity 
equations cannot be obtained from the chemical kinetics. In this case, the production 
terms in the energy equation are eliminated by formulating this equation in terms of 
enthalpy. The energy equation in terms of total enthalpy is 


p ( u 11 + v \ - v 15. + pKU v 
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( 8 ) 
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where 


H- 

2 


(9) 


The production terms that appear in the species continuity equations are eliminated 
by introducing the concept of elemental mass fractions. As long as no nuclear reactions 
occur, the elemental mass fractions remain fixed and unchanged during chemical reactions. 
The relation between the elemental and species mass fractions (the Shvab- Zeldovich trans- 
formation (ref. 22)) is given by 


/V 




8 u 


Mf 


( 10 ) 


The elemental continuity equations for the elements can be obtained by multiplying equa- 
M*C. 

tion (5) by 5. — 8 — 1 and summing over i. The resulting elemental continuity equations 

are 
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'V 


'Xy . 

bC 


_ + v_ 
^1 + n k bs an 




(1 + n/c) (r + n cos 6 )J Un 


[(1 + n/<)(r + n cos 6 $ J g ] 


(ID 


where 


K 



5. 

n 


M* 

— J i 
M| 1 


( 12 ) 


Introducing the elemental mass fraction not only eliminates the production terms in equa- 
tions (4) and (5) but also reduces the number of equations to be solved, since there is one 
equation (eq. (11)) for each element rather than one equation (eq. (5)) for each species. 

Equations (1) to (7) are the governing relations for frozen and nonequilibrium flow. 
Equations (4) and (5) are replaced with equations (8) and (11), respectively, for equilibrium 
flow. In these equations, additional relations are necessary to specify mixture quantities 
and the diffusion mass fluxes. Appendix B presents the relations that are used to calculate 
the thermodynamic and transport properties for each chemical species and the mixture 
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transport properties. The values for the remaining mixture quantities are determined by 
the following relations: 

Frozen specific heat: 


N 


s ■ L c i c p,i 

i=l 


where 


(13) 


c = p 1 

i p 


(14) 


Molecular weight: 


M*=! 


N s 

Z- 

w M ! 


(15) 


Enthalpy: 
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h * E h i c i 

i=l 


(16) 


where h. includes the enthalpy of formation of the ith species. 

The mass flux due to concentration gradients can be written as (ref. 1) 
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where 
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( 18 ) 
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and 


NLe,i 
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L-y Mf 
3 = 1 1 
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j = 1 1 3 3 
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(19) 


In equations (18) and (19), are the multicomponent Lewis numbers and are the 

binary Lewis numbers. The relative mass flux for the elements can be written as 
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where 
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Ab ik = 
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( 22 ) 


and L is an arbitrary constant (generally taken as one). Equations (17) and (20) can be 
simplified to 


and 


M 






(23) 


(24) 


respectively, for binary diffusion. Note that K is equal to either ? , a constant binary Lewis 
number, or J.j, a variable binary Lewis number. 
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Boundary Conditions 

Conditions at the body surface .- The no- slip boundary conditions are used in this 
study. The surface conditions for rj = 0 (see fig. 2) are 


u - 0 



where the mass-injection rate 


(25) 

(26) 


m = (pv) 


(27) 


is specified. In the calculations presented in this paper, the mass-injection distribution 
along the body surface is assumed to be 


m = m 0 e 


-3.47r2 


(28) 


where m 0 is a specified stagnation- point mass-injection rate. This distribution is used to 
demonstrate the effect of mass injection on flow-field properties and to demonstrate the 
capability of the present analyses to account for distributed mass injection. The same 
distribution is used in reference 8. 

The wall temperature for this study is specified as 


T . = Constant 
w 


(29) 


The surface total enthalpy is given as 
N s 

H = h i C i (30) 

i^l 

The surface boundary conditions for the elemental and species equations are derived as 
follows. (See fig. 2.) The species i are transported away from the surface at the rate J. 
by diffusion and at the rate pvC i by convection. At the same time, the species i are being 
convected to the surface at the rate pvCj_. This flux may be considered to be that from 
an ablating surface or the injectant from a porous solid. The surface species concentra- 
tions are given by 

pvC.+ e 2 J. = ( pvCj)_ (31) 
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Since 


pv = (pv) p m 


(32) 


equation (31) may be written as 

(C.- C._)m= -e 2 J. (33) 

and in terms of elemental boundary conditions as 

ft -eym = -, 2 J ( (34) 

Conditions at the shock .- The conditions imposed at the shock are calculated by 
using the Rankine-Hugoniot relations. The flow is assumed to be either in chemical equi- 
librium or frozen at the free -stream composition. The nondimensional shock relations 
are as follows: 

Mass: 

P s v s = -sin a (35) 


u^ = cos a 

Momentum: 


Energy: 


1 .2 ft 1 

p = + sin a [ 1 - _ 

7 TC M 2 \ P s 


u 1 sin 2 a 
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P T R* 
= ^ S 


s M*C* 
s p,® 


Enthalpy: 


N. 


h s h i,s^i,£ 


i= 1 


(36) 


(37) 


(38) 


(39) 


(40) 
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Elemental composition: 


%,s = V (41) 

and when the chemistry is frozen across the shock, 

Species composition: 



(42) 


Surface transport .- The rate at which heat q is transferred to the surface boundary 
is given by the sum of the convective, conductive, and diffusion contributions. The net 
nondimensional heat transfer to the solid interior is given by 


q 


/ N s 

M.kK* t 

\ ' •» k 



N 





J w 


The Stanton number is given by 

q 

N St = " h -H 

n oo n W 

and the skin-friction coefficient is given by 


(43) 


(44) 


C. = 2?(n**S (45) 

\ 3n /w 

A transformation is applied to the previous nondimensional viscous shock- layer 
equations and boundary conditions to simplify the numerical computations. The transfor- 
mation relations and the transformed equations and boundary conditions are given in 
appendix C. 


Chemical Composition 

Analyses of chemically reacting flows are frequently simplified by considering the 
limiting cases of frozen and equilibrium chemical behavior. However, the question 
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naturally arises: Which of these two limiting cases is more descriptive of reality? Of 
course, to answer this question the detailed finite-rate behavior of the chemical reactions 
must be accounted for. An analysis of a complex chemically reacting gas mixture, where 
the chemistry is nonequilibrium, presents problems that are not associated with either 
frozen or equilibrium treatments. These problems are primarily those of defining a 
realistic reactions model, of obtaining the appropriate rate constants, and of obtaining a 
numerical solution. In general, it is the latter problem that precludes a nonequilibrium 
treatment of reacting flows. 

In this study, the chemical reactions will be confined to a system of neutral air 
species (O, Og, N, Ng, and NO). However, the chemical models are not limited to this 
system of species. 

Equilibrium .- The equilibrium chemical composition is determined (for a given 
temperature, pressure, and elemental composition) by using a free- energy minimization 
analysis that is described in reference 23. Since the criterion for equilibrium at constant 
temperature and pressure is that the change in free energy be zero, dF* = 0, the equilib- 
rium composition is determined when the total free energy of the mixture is made a 
minimum with respect to any possible change in composition. Details concerning the 
equilibrium calculation can be found in reference 23. 

Nonequilibrium.- When chemical reactions proceed at a finite rate, the rate of 

production terms w^ are required. The production terms appear in the energy equation 

(eq. (4)) when formulated in terms of temperature and in the species continuity equations 

(eq. (5)). For a multicomponent gas with N„ reacting chemical species and N chemical 

s r 

reactions, the chemical equation describing the overall change from reactants to products 
may be written in the general form 



v x t 





(46) 


where r = 1, 2, . . ., N r and Nj is equal to the sum of the reacting species (Ns) plus the 
number of catalytic bodies. The quantities a. and 0. are the stoichiometric coefficients 
for reactants and products, respectively, whereas and k* r are the forward and back- 
ward rate constants. The quantities X? denote the concentrations in moles per volume. 

The catalytic bodies (Nj - N s ) may be chemical species or linear combinations of species 
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that do not undergo a chemical change during the reaction. The rate of change of any 
species as a result of a particular reaction (ref. 24) is 






(47) 


Equation (47) may be rewritten as 


(w*) r 


M* 




v 


) 



3=1 



(48) 


where the mole-mass ratio 7? is defined as 


and 


y^< 


X * c. 

3 = 3 

p * M| 


N 




i=l 


dX* dp* 
w 5 ? - 1 - 1 


1 1 dt* dt* 


(3 = 1, 2, . . ., N s ) 




F. Z (]-N s ),i T i G “ N s+1> • • •> N j) 


(49) 


(50) 


The constants j are determined from linear dependence of the catalytic bodies 

upon the N species. Values of these constants are given in reference 11 and in table I. 
s 

In order to find the net mass rate of production of the ith species per unit volume, 
equation (48) must be summed over all reactions r. Thus, one obtains 


N r 

N i 


N j 


*t = M i L - V> 


” ( 7 ^P*) aj,r - k* r 


('y ^ p ^ ^ 

r=l 

„ 3=3 


==•] 

— 


which is the general rate equation for a gas mixture. 


(51) 
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The chemical reactions used in the study are as follows: 


r = 1 

0 2 + — 20 + 

r = 2 

n 2 + m 2 ^:2n+ m 2 

r = 3 

N 2 +N^2N+N 

11 

NO + Mg ^ N + O + 

r = 5 

NO + O 0 2 + N 

r = 6 

N 2 + O NO + N • 


y 


( 52 ) 


The reaction constants for these equations are expressed in the modified Arrhenius form, 
where the forward rate is given as 


C2* ( Cl* x 10° \ , / . \ “ a r 

jj* T* r exp log e CO — — — ' 1 i mole \ 




T* 


s\ 3 
\cm 


(53) 


and the backward rate is given as 

D2 r 

k b,r = T * exp | log e DO r 


where 


and 


D1 ? x 10 ) 1 1 

T* Mcm ! i 


N i 

“r= E V - 1 

i=l 


n j 

«r= E V - 1 
i=l 


(54) 


(55) 


(56) 


Values for the coefficients in equations (53) and (54) are taken from the compilation of 
experimentally determined rate constants given in reference 25 and are tabulated in 
table II. Therefore, for a specified temperature, density, and species composition, equa- 
tions (51) to (54) are used to determine the production rate of a multicomponent gas. The 
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manner in which the production terms are written for numerical solution is discussed in 
appendix C. (See eqs. (C17) and (C19).) 

Method of Solution 

The procedure for solving the viscous-shock-layer equations is presented herein. 
First, the finite- difference expressions used to transform the differential expressions 
to algebraic expressions are presented. Then the solution procedure is discussed. 

Finite- difference expressions .- The derivatives are converted to finite-difference 
form by using Taylor's series expansions. A variable grid spacing (fig. 3) is used in the 
r?-direction so that the grid spacing can be made small in the region of large gradients. 
Three-point differences are used in the 17 -direction, and two-point fully implicit differences 
are used in the f -direction. Truncation terms of order a f (first order accuracy) and 
either A *? n ^ or ( A *7 ~ A 7 , n _i^ (second order accuracy) are neglected. Atypical 
finite -difference expansion of the standard differential equation (see appendix C, 
eq. (C 8 )) gives 


A W. 


B W 


n m,n-l + n " m,n 


C W m - = D 
n m,n+l n 


(57) 


The coefficients A, B, C, and D are used to represent the coefficients after the finite - 
difference expansion of equation (C 8 ). The subscript n denotes the grid points along a line 
normal to the body surface, whereas the subscript m denotes the grid stations along the 
body surface. Equations (57) along with the boundary conditions constitute a system of 
the tridiagonal form, for which efficient computational procedures are available. (See 
ref. 26.) 

Overall solution procedure .- For specified free-stream conditions and body geometry, 
a stagnation streamline solution is obtained. With the stagnation streamline solution pro- 
viding the initial conditions, the conditions at the shock providing the outer boundary con- 
ditions, and the conditions at the wall taken as the inner boundary conditions, the numerical 
solution is marched downstream to the desired body location £. The first solution pass 
provides only an approximate flow-field solution, because the following assumptions are 
used in the first solution pass: 

(a) The thin- shock- layer form of the n-momentum equation (eq. (C15b)) is used 

(b) The stagnation streamline solution is independent of downstream influence 

(approximation of local similarity where n 2 s = 0 ) 

(c) The term dn s /d£ is equated to zero at each body station 

(d) The shock angle a is assumed to be the same as the body angle 0 . 


20 



These assumptions are then removed by making one or more additional solution 
passes. For the current study, a total of two solution passes are used since the two passes 
resulted in a converged flow-field solution. For the second solution pass, the thin- shock- 
layer form of the normal momentum equation (eq. (Cl 5b) is replaced with equation (Cl5a).) 
The v component of velocity that is used in equation (Cl 5a) is the average of the first and 
second pass values. Also, once the first solution pass has been computed, the values of 
n 2s and dn s /d| are calculated and used in the second solution pass to remove approxima- 
tions b, c, and d. Hence, the viscous- shock- layer equations axe solved as parabolic equa- 
tions, and yet retain effects which are elliptic and hyperbolic in nature. This solution 
procedure is programed for the Control Data 6600 computer. 

Shock solution. - The shock solution procedure at any location is identical for the 
first and subsequent solution passes. However, the shock angle a is defined differently 
for the first and subsequent solution passes. The shock angle a is set equal to the local 
body angle 0 for the first solution pass. For subsequent solution passes, the shock angle 
is defined as 


i n L 

a = 6 + tan" 1 — (58) 

l4 “, 

Figures 4(a) and 4(b) present the solution flow diagrams for equilibrium and frozen shocks, 
respectively. 

Solution procedure at station m .- The viscous- shock-layer equations are solved at 
any body station m (see fig. 3) in the order shown in figure 5(a) for equilibrium flow and 
in the order shown in figure 5(b) for either frozen or nonequilibrium flow. The governing 
equations are uncoupled and the dependent variables are solved one at a time in the order 
shown in figure 5. First, the shock conditions are calculated to establish the outer boundary 
conditions. Then the converged profiles at station m - 1 are used as the initial guess for 
the profiles at station m. The solution is then iterated locally until convergence is achieved. 
For the stagnation streamline (m = 1), guess values for the profiles are used to start the 
solution. 

Each of the second- order partial differential equations are individually integrated 

numerically by using the tridiagonal formalism (eq. (57)). The global continuity equation 

is used to obtain both the shock standoff distance and the v components of velocity. By 

integrating equation (C14) between the limits of n = 0 and n = 1 at station m, an implicit 

equation for n is obtained. For the v component of velocity at i?, equation (C14) is inte- 
s 

grated with respect to n between the limits of 0 to n . The pressure p is determined at 
station m by integrating the normal momentum equations (Cl 5) with respect to i? between 
the limits of 1 to r? . The equation of state is used to determine the density. 
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VERIFICATION OF ANALYSIS 


In order to assess the overall results obtained for the present analyses, compari- 
sons are made of the results for each of the three chemistry models with data reported in 
the literature. Table HI gives a summary of the viscous-shock-layer and boundary- layer 
analyses for which comparisons are made with the present analyses. 

All computations employed a constant step size of 0.2 in the ^-direction and a 
variable step size in the v -direction. A total of 50 grid points is used at each body sta- 
tion m. Also, two solution passes are made in the £ -direction. The two solution passes 
provide a converged shock shape and flow-field solution for the problems considered in 
this study. Table IV summarizes the altitude, velocity, and free-stream properties used. 
The free-stream properties are those given by the U.S. Standard Atmosphere (see ref. 27). 

For each body station, the convergence criterion is that the relative difference 
between the current and previous iteration values be less than 0.001 for both the tempera- 
ture and tangential velocity derivatives at the wall. This criterion is used for all calcula- 
tions except for those where the temperature gradient approaches zero because of mass 
injection. For this condition, only the velocity derivative is checked. A minimum of three 
interations is required at each body station. 

Frozen Flow Comparisons 

The frozen flow calculations are compared with the data of Whitehead (ref. 8). Even 
though the solution procedure used herein is essentially the same as that used by White- 
head, several differences exist. Whitehead's results account for shock slip at the outer 
boundary condition and slip and temperature jump at the surface, whereas the present 
analysis uses the Rankine-Hugoniot relations to describe the shock conditions and assumes 
no slip or temperature jump at the wall. Furthermore, the expressions used to calculate 
the thermodynamic and transport properties are not the same as those used by Whitehead. 

Figure 6 shows comparisons of the present analysis with that of Whitehead's for 
shock standoff distance, skin-friction coefficient, and nondimensional heat-transfer rate 
when the chemistry is frozen at the free-stream composition. The heat-transfer results 
(fig. 6(c)) are in excellent agreement; however, the agreement is not as good for shock 
standoff distances (fig. 6(a)) or skin-friction coefficients (fig. 6(b)). Table V presents a 
detailed listing of the shock and wall values that result from the present calculation. 
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Equilibrium Flow Comparisons 


The equilibrium results are compared with the results obtained from six previous 
studies. Figure 7 compares the present results with those obtained by Whitehead (ref. 8). 
The Lewis number for comparison purposes is assumed to be one. Whitehead's equilibrium 
analysis assumes air to be one "effective" species with a constant molecular weight equal 
to the free- stream value, whereas the present analysis treats air as a multicomponent 
mixture of five species. The present calculation predicts greater values for shock stand- 
off distance (fig. 7(a)) (as would be expected because of Whitehead's assumption concerning 
the molecular weight), skin-friction coefficient (fig. 7(b)), and heat transfer (fig. 7(c)). 

Table VI(a) presents a detailed listing of the stagnation shock and wall values for the 
present calculation. 

Results of the equilibrium air calculations are also compared with the stagnation 
results of Edelman and Hoffman. (See ref. 15.) Comparisons of stagnation temperature 
and species concentration profiles are presented in figure 8. Figures 8(a) and 8(b) show 
a comparison of the temperature and species profiles, respectively. Table VI(b) presents 
a detailed listing of the stagnation shock and wall values for the present calculation. 

Figures 8(c) and 8(d) show the same comparisons for a different set of flow conditions. 
Table VI(c) presents the stagnation shock and wall values for the present calculation. 

When the results of the two calculations for both sets of conditions are compared, 
two differences are apparent: first, the shape of the temperature profiles, and second, 
the shock-layer thickness. The curvature of the temperature profiles as calculated by 
reference 15 decreases monotonically with increasing distance from the surface. Tem- 
perature profiles resulting from the present equilibrium calculations are not smooth, 
especially at the higher altitude (lower pressure) conditions. However, it is believed that 
the temperature profiles resulting from the present calculations can be substantiated by 
considering the chemical reactions that occur within the shock layer. A more detailed 
discussion of this assessment is given in the next section. 

The differences in shock-layer thickness are expected since the results presented 
in reference 15 were for the stagnation streamline only and did not account for the down- 
stream influence on shock standoff distance. When the downstream influence on shock 
shape is accounted for as in the present analysis, a larger shock- layer thickness results. 
This effect was shown in reference 6 for an ideal gas. 

Figure 9 shows comparisons of the present equilibrium analysis with the stagnation 
viscous-shock-layer analysis of Goldberg and Scala (ref. 16). Comparisons of static 
enthalpy, normal velocity, and tangential velocity are presented. The calculations are for 
a shock Reynolds number of 100. V/ith the exception of shock standoff distance, the 
results of the two analyses are in good agreement, particularly near the wall where each 
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of the three profiles are in excellent agreement. Table VT(d) gives the shock and wall 
conditions for the calculation. 

As a final comparison of the equilibrium calculations, the viscous- shock- layer 
results are compared with boundary-layer solutions of Anderson and Lewis, Blottner, and 
Smith as presented in reference 28. The calculations are for flow over a 10° half-angle 
hyperboloid, A detailed listing of stagnation shock and wall conditions are given in 
table VT(e). 

The pressure distribution for the boundary- layer calculations was obtained from 
modified Newtonian theory. This distribution is in excellent agreement with the pressure 
distribution calculated with the viscous- shock- layer analysis. Therefore, comparison of 
the two pressure distributions is not presented. 

Figure 10 shows comparisons of skin-friction distributions. The three boundary- 
layer solutions and the viscous- shock- layer solution are in good agreement. The agree- 
ment for Stanton number distributions is even better, and therefore, comparisons are not 
presented. 

The aforementioned comparisons of the present equilibrium viscous shock- layer 
solutions with previous shock- layer and boundary -layer solutions have generally shown 
good agreement. Where apparent differences exist, such as shock standoff distance, the 
differences can be rationalized. Note that the comparisons encompass a wide range of 
flow conditions with shock Reynolds numbers ranging from 100 to 17 931. 

Nonequilibrium Comparisons 

Results of the nonequilibrium calculations are compared with the stagnation viscous- 
shock- layer results of Blottner. (See ref. 11.) Comparisons are made of temperature, 
tangential velocity, and species concentration profiles. 

For the chemical kinetics model, Blottner considered two additional species (NO + 
and e") and one additional reaction equation (N + OZT NO + + e"); however, for the remain- 
ing reactions, identical rate expressions are used. The calculations are for an equilibrium 
catalytic wall. Both analyses use the Rankine -Hugoniot relations to describe the shock 
conditions and use the no- slip or temperature jump conditions at the wall. Also, for pur- 
poses of comparison, the present calculation is made with a constant Lewis number of 1.4. 

For two reasons, differences in shock- layer thickness are expected. First, the 
results of reference 11 do not account for the downstream influence of the stagnation 
shock standoff distance as does the present analysis. When the downstream effect is 
neglected, a significant underprediction of stagnation shock standoff distance can occur. 
(See refs. 6 and 29.) Second, the analysis of reference 11 did not include the effect of 
shock-layer thickness on the shock radius of curvature as did the present analysis. 


24 



Kaiser and Flugge-Lotz (ref. 29, p. 91) found that when the shock-layer thickness con- 
tribution to the radius of curvature is neglected, an underprediction of the stagnation 
shock standoff distance occurs. 

Table VII presents a detailed listing of stagnation shock and wall conditions for the 
present calculation (a constant Lewis number of 1.4). Figures 11(a) and 11(b) show com- 
parisons of temperature and velocity ratio profiles, respectively. The difference in shock 
standoff distance is readily apparent. These results also show that the temperature and 
velocity ratio profiles are in good agreement near the wall. 

Figure 11(c) shows a comparison of the species concentration profiles for the two 
analyses. Appreciable differences exist. These differences are due primarily to differ- 
ences in shock standoff distance. The shock standoff distance for the present analysis 
is approximately 22 percent greater than that of reference 11. Consequently, a larger 
distance is available for reactions to occur - in particular, the dissociation reactions. 
This situation is evident in figure 11(c) which shows that for a given value of n, the amount 
of dissociated species predicted by the present calculation is greater than that predicted 
by reference 11. 


Significance of Comparisons 

The comparisons of the present solutions with previously reported solutions 
(table m) for frozen, equilibrium, and nonequilibrium flows show generally good agree- 
ment. Yet none of the analyses in table III are as detailed as the present analyses. 
Therefore, for problems similar to those that have been discussed, the more approximate 
solutions will provide good accuracy. However, this conclusion will not be true for many 
flow conditions. For example, the importance of including multicomponent diffusion 
becomes more significant when the species molecular weight disparity becomes greater 
than that of the neutral air species. (See ref. 4.) Furthermore, it is important to note 
that the present analyses are more general than those listed in table III. The present 
analyses provide stagnation and downstream solutions, account for multicomponent diffu- 
sion, account for mass injection, and provide solutions where the flow-field chemical 
species are not limited to the neutral air species. 

RESULTS AND DISCUSSION 

Numerical solutions to the viscous -shock- layer equations are presented and dis- 
cussed herein. The solutions are those for flow about a hyperboloid (total included angle 
of 45°) at an angle of attack of 0°. Results are presented which demonstrate the effects 
of chemistry model, diffusion model, surface catalycity, and mass injection on flow 
parameters and surface transport. 
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All computations are made with the numerical step-size values and convergence cri- 
teria that were presented in the previous section. Also, all the subsequent results will be 
for a free-stream velocity of 6.10 km/s, an altitude of 60.96 km in the earth's atmosphere, 
a 2.54-cm nose radius, and a surface temperature of 1500 K. These flow conditions are 
representative of those experienced by slender vehicles and the leading edges of manned 
entry vehicles when entering the atmosphere from a low earth orbit. Note, however, that 
the analyses are not restricted to the slender bodies considered in the present study. 

Chemistry Models 

In this section results of calculations using frozen, equilibrium, and nonequilibrium 
chemistry are presented. Multicomponent diffusion, no mass injection, and a noncatalytic 
wall are assumed. First, velocity, temperature, and species concentration profiles are 
presented for each chemistry model. Next, the equilibrium temperature profiles are shown 
and their dependence on chemical reactions and pressure variations are discussed. Then 
the effect of each chemistry model on tangential velocity, temperature, and chemical 
species profiles, heat transfer, Stanton number, skin friction, shock standoff distance, 
and wall pressure is shown. 

Figures 12 to 15 show the velocity, temperature, and species profiles at various 
body stations for each chemistry model. There are several features of these profiles 
that are of interest. The tangential velocity profiles (figs. 12(a) to 12(c)) are practically 
linear in the outer region of the shock layer; that is, the flow has an outer region of essen- 
tially uniform vorticity. This result is in marked contrast to the classical velocity 
boundary- layer profile that is derived with a zero velocity gradient imposed at the outer 
boundary. The flow considered here is in the viscous-layer regime (ref. 30) where the 
viscous effects extend over an appreciable fraction of the shock layer and thereby invali- 
date the boundary- layer concept. 

The equilibrium temperature profiles are shown in figure 13(b). Near the wall, the 
equilibrium temperature profiles do not show the monotonically decreasing curvature 
associated with nonequilibrium (fig. 13(a)) and frozen (fig. 13(c)) profiles. The equilibrium 
temperature profiles have two distinct inflection points. From the body to the shock, the 
profiles are characterized by concave, convex, and concave segments. 

The bulges appearing in the temperature profiles can be explained in terms of 
chemical reactions. This effect is demonstrated in figure 16, where stagnation tempera- 
ture profiles and species concentration profiles are shown. Temperature profiles for an 
equilibrium chemistry solution and a frozen chemistry solution (where the chemistry is 
frozen at the equilibrium shock composition) are shown. The temperature profile for the 
equilibrium calculation is much fuller than the profile for the frozen calculation. This 
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condition should exist, of course, because exothermic reactions occur near the wall. 
Essentially all chemical reactions occur within that half of the shock layer which is adja- 
cent to the wall; that is, essentially all the atom recombinations occur near the wall. The 
recombinations also occur over two distinct and separate regions. The oxygen recombines 
between the wall and r? = 0.11. This condition causes the inner bulge in the temperature 
profile. Then there is a small region (rj of 0.11 to 0.16) where the recombination process 
is negligible. Finally, practically all the atomic nitrogen recombines between v values of 
0.16 and 0.5. This recombination causes the outer bulge in the temperature profile. 

The extent to which the equilibrium temperature profiles differ from smooth profiles 
will depend upon the pressure, and hence the altitude for a given free-stream velocity. 

This effect is demonstrated by Hansen in figure 1 of reference 31. The data of Hansen 
show that the recombination reactions occur over narrow temperature intervals at low 
pressures (0.0001 to 0.1 atm); however, the temperature intervals required for nitrogen 
and oxygen recombination increase significantly with increasing pressure. Furthermore, 
the temperature interval between the completion of nitrogen recombination and the onset 
of oxygen recombination decreases with increasing pressure. Consequently, as the pres- 
sure increases, the chemical composition experiences a more gradual change with tem- 
perature. The pressures are 0.1 atmosphere or less for the temperature profiles shown 
in figure 13(b); therefore, recombinations over narrow separated temperature intervals 
would be expected. 

Figure 14 shows the nonequilibrium species profiles at different body stations for a 
noncatalytic wall. The chemistry across the shock is frozen at the free-stream composi- 
tion (Cq = 0.24; C ]$2 = 0.76). The figures show that most of the dissociation occurs in the 
outermost part of the shock layer. Furthermore, only a small amount of recombination 
occurs near the wall. This effect is due to the low-density flow conditions which reduce 
the effectiveness of recombination reactions; consequently, the amount of dissociated flow 
at the wall is significant (approximately one-fifth of the total mass). This result is 
important because the energy invested in dissociation is not recovered. 

Figure 15 shows the equilibrium concentration profiles at different body stations. It 
is seen that a large percentage of flow is dissociated in the outer part of the shock layer. 
The amount of dissociation decreases rapidly downstream. Essentially all the dissociated 
species recombine in a relatively small region near the wall, with the wall composition 
being approximately that of the free-stream composition. This statement, of course, 
means that the energy of dissociation is recovered. The significance of this statement 
in terms of heat-transfer rate will be discussed. 

Comparisons of results using the three chemistry models are shown in figures 17 
to 23. The effects on tangential velocity and temperature profiles, shock standoff distance, 
heat-transfer rate, Stanton number, skin friction, and wall pressure are shown. Note that 
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all these quantities are significantly influenced by the chemistry model with the exception 
of the tangential velocity profiles (fig. 17) and the wall pressure distribution (fig. 23). Of 
particular significance are the large differences in shock- layer temperatures for the dif- 
ferent chemistry models. Figures 18(a) and 18(b) show the temperature profiles at an £ of 
0.0 and 2.0, respectively, for each chemistry model. The conditions at the shock were 
identical for the frozen and nonequilibrium calculation since the chemical composition was 
assumed to be frozen across the shock at the free- stream composition. For the stagnation 
streamline (fig. 18(a)), the frozen shock temperature is more than twice that of the 
equilibrium shock temperature. Moreover, the equilibrium temperatures across the shock 
layer are much less than those of the frozen or nonequilibrium values except in the region 
near the wall. This condition is due to the large amount of energy invested in dissociation 
rather than in static temperature. Also, the nonequilibrium temperatures are less than the 
frozen values because some dissociation occurs within the nonequilibrium shock layer. At 
the downstream station, £ = 2.0 (fig. 18(b)), the temperature differences are noticeably less 
than those at the stagnation station because the equilibrium and nonequilibrium composi- 
tions are more nearly like the frozen composition at the lower temperatures. 

Figure 19 shows how the chemistry model influences shock standoff distance. The 
shock standoff distances corresponding to the equilibrium calculation are considerably 
less than those for frozen and nonequilibrium calculations. Also, the shock standoff dis- 
tances for the nonequilibrium calculations are somewhat less than those for the frozen 
calculations. These results are expected since the shock- layer densities for the non- 
equilibrium and especially the equilibrium calculations are greater than the frozen flow 
densities. 

Consider next the effect that the chemistry model has on heat-transfer rate distribu- 
tion. As is shown in figure 20, the chemistry model has a very pronounced effect on non- 
dimensional heat transfer. Both the frozen and equilibrium chemistry models significantly 
overpredict the heat transfer to a noncatalytic wall. This condition is especially true in 
the stagnation region where the frozen and equilibrium values are approximately 50 per- 
cent greater than the nonequilibrium values. This fact is particularly significant because 
the entry corridor for slender entry vehicles includes the altitude -velocity conditions used 
in the present calculations. Furthermore, the heating rates are significant at these condi- 
tions, and therefore, the need for nonequilibrium analyses is clearly evident from these 
results. 

When the heat-transfer rate is expressed in terms of Stanton number (fig. 21), the 
results for the three chemistry models are in closer agreement than were the non- 
dimensional heat-transfer rates (fig. 20). As was pointed out previously (fig. 14), a sig- 
nificant amount of energy is invested in dissociation at the wall in nonequilibrium flow. 
Consequently, the enthalpy at the wall is larger for the nonequilibrium and noncatalytic 
wall calculation and the enthalpy potential is smaller. The maximum difference between 



the equilibrium and nonequilibrium Stanton numbers is about 13 percent and occurs at the 
stagnation point. 

The effect of the chemistry model on skin-friction distribution is shown in figure 22. 
This figure shows that the skin- friction distributions for the equilibrium and nonequilibrium 
calculations are about the same. The values for the frozen calculation are greater than the 
nonequilibrium values, but the differences never exceed 9 percent. 

The effect of the chemistry model on wall pressure distribution is shown in figure 23. 
Also shown is the Newtonian pressure distribution. The nonequilibrium and frozen 
pressure distributions are identical. The equilibrium and Newtonian pressure values are 
in close agreement and are less than the nonequilibrium pressure. At £ = 3.0, the non- 
equilibrium pressure is 10.8 percent greater than the equilibrium and the Newtonian 
pressures. 


Catalytic Wall Effects 

Results presented in the previous section demonstrated how the flow-field chemistry 
can influence flow parameters and surface transport. This section demonstrates how sur- 
face catalysis influences the same flow parameters and surface transport. Either non- 
catalytic or equilibrium catalytic walls are assumed. For the equilibrium catalytic wall, 
the gas composition at the wall is the equilibrium composition for the wall temperature 
and pressure. The elemental composition at the wall is governed by multicomponent dif- 
fusion and is, therefore, not necessarily equal to the free- stream elemental composition. 

Comparison of the noncatalytic and equilibrium catalytic wall calculations is made. 
The diffusion for both calculations is multicomponent. Table VII provides a detailed list- 
ing of the stagnation shock and wall values for the noncatalytic and equilibrium catalytic 
calculations. 

The wall catalytic activity has negligible effect on temperature and velocity profiles, 
and therefore comparisons are not presented. However, the species profiles for the flow 
configuration examined are necessarily influenced by the wall catalytic activity. This 
influence is shown in figures 24(a), 24(b), and 24(c), where comparisons are made at an £ 
of 0, 1.0, and 3.0, respectively, of nonequilibrium species profiles with noncatalytic and 
equilibrium catalytic wall conditions. It is seen that the species profiles for the inner- 
most 70 percent of the shock layer are influenced by the wall catalytic activity. 

Figures 25 to 28 show how wall catalytic activity influences shock standoff distance, 
heat transfer, Stanton number, and skin friction. With the exception of heat transfer, the 
wall catalytic activity has little effect. The maximum relative differences between the 
catalytic and noncatalytic calculations are 3.5 percent for shock standoff distance, 

4.9 percent for skin friction, and 9.2 percent for Stanton number. The maximum relative 
difference in heat transfer is 48 percent. 
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For purposes of comparison, figure 26 also includes the heat-transfer distribution 
for equilibrium flow. It is seen that the results for nonequilibrium flow with an equilibrium 
catalytic wall and the results for equilibrium flow are in close agreement. For both of these 
calculations, the energy invested in dissociation is recovered. However, the results for 
nonequilibrium flow with a noncatalytic wall show an appreciable amount of dissociation 
at the wall (see fig. 14), and, hence, the wall heat transfer is substantially reduced. 

Diffusion Models 

Prior to this study, stagnation and downstream viscous -shock-layer results that 
treated the diffusion as multicomponent were not available. In this section, comparisons 
are shown between multicomponent and binary diffusion models. The effects on flow 
parameters and surface transport are- presented for both equilibrium and nonequilibrium 
air. 

One of the binary approximations is that all species have the same diffusion coeffi- 
cient that is equal to the coefficient for molecular nitrogen diffusing into atomic oxygen. 

This approach is logical because atomic oxygen and molecular nitrogen are the dominant 
species in the shock layer. (See figs. 14 and 15.) The multicomponent diffusion coeffi- 
cients are determined according to equation (B15). Note that this approach for evaluating 
the binary diffusion coefficient results in a variable Lewis number. 

No mass injection and a noncatalytic wall are assumed. With the exception of the 
species profiles, essentially no effect of diffusion model is observed. Even the effect on 
species profiles is small as is shown in figures 29 and 30. 

Figures 29(a) and 29(b) show comparisons of binary and multicomponent species 
profiles for equilibrium air at an £ of 0.0 and 1.0, respectively. More molecular oxygen 
and less molecular nitrogen are present at the wall for multicomponent diffusion. 

The elemental composition is not a constant for multicomponent diffusion as is the 
case for binary diffusion. However, the departure from elemental invariance is small for 
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the equilibrium air calculation. The elemental composition at the wall is C Q = 0.278 and 
C N = 0.722 at £ = 0.0, and C Q = 0.279 and C N = 0.721 at £ = 1.0, as opposed to the free- 
stream values of Cq = 0.24 and *= 0.76. 

Figures 30(a) and 30(b) show comparisons of species profiles for multicomponent 
and binary diffusion models at £ values of 0 and 1.0, respectively. The flow is non- 
equilibrium air. Once again the diffusion model has a small effect on species profiles. 

The elemental wall concentrations are Cq = 0.221 and 6^ = 0.778 at £ =0, and Cq = 0.222 
and C N = 0.779 at £ = 1.0. 

The other flow parameters and surface transport values are essentially the same 
for both diffusion models, and therefore comparisons of results are not presented. For 
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example, the differences in stagnation heat-transfer rates, between binary and multi- 
component values, were 0.74 percent for nonequilibrium flow and -0.84 percent for 
equilibrium flow. 

The question of whether a variable binary Lewis number, as previously described, 
provides better results than a constant Lewis number is now considered. For air, the 
most frequently used values of Lewis number are 1.0 and 1.4. Calculations are made for 
both equilibrium and nonequilibrium air (noncatalytic wall) with constant Lewis numbers of 
1.0 and 1.4. For the equilibrium calculations, a Lewis number of 1.0 underpredicts and a 
Lewis number of 1.4 overpredicts the heat-transfer rate. The differences in stagnation 
heat-transfer rates were 10.3 percent for a Lewis number of 1.4 and -5.5 percent for a 
Lewis number of 1.0. For the nonequilibrium calculations, the effects of diffusion model 
on heating rates were negligible since the errors were less than 0.8 percent. Therefore, 
these results show that the diffusion model has a small effect on most flow-field proper- 
ties and surface-transport values. However, differences as large as 10 percent can be 
incurred for equilibrium heat-transfer rates when a constant Lewis number of 1.0 or 1.4 
is used in lieu of multicomponent diffusion. Also, the differences introduced by using a 
binary diffusion approximation can be significantly reduced by using a variable Lewis 
number where the diffusion coefficient is based on the two dominant species. 

Mass Injection 

In this section, the results of computations with mass injection are considered. In 
all cases the mass injection rate distribution is given by equation (28). Figure 31 shows 
the injection distribution for different values of the nondimensional stagnation injection 
parameter m 0 . It is seen that the injection rate approaches zero at a £ of about 1.5. 

Figures 32 to 39 and table VIII show how injecting equilibrium air into nonequilibrium 
air influences flow parameters and surface transport. Multicomponent diffusion and a 
noncatalytic wall are assumed. Stagnation mass injection rates as large as 0.4 are 
considered. 

Figure 32 shows comparisons of species profiles for no injection and an injection 
rate of 0.2. These comparisons are made at a | of 0, 1.0, and 3.0. At the stagnation 
point, the shock- layer chemical composition is altered substantially because of mass 
injection (fig. 32(a)). At the wall, the mass of dissociated flow is 6 percent for a 0.2 
injection rate and 24 percent for no injection. Downstream the effect of mass injection 
on chemical composition decreases (fig. 32(b)) because the injection rate is smaller 
(eq. (28)). At £ = 3.0, which is beyond the body station where mass injection is zero, the 
chemical composition of the two shock layers (fig. 32(c)) is about the same. 

Figures 33 and 34 show the stagnation tangential velocity and temperature profiles, 
respectively. These results qualitatively are as one would expect. That is, the shock- layer 
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thickness increases and the tangential velocity and temperature gradients decrease with 
increasing mass injection rate. 

Figure 35 shows the effect of mass injection on shock standoff distance. For stag- 
nation injection rates of 0.2 or greater, the injectant increases the shock standoff distance 
as far downstream as the calculations are made ( £ = 3.0). 

Figure 36 shows the effect of injecting equilibrium air into nonequilibrium air on 
nondimensional heat-transfer distribution. As the injection rate increases, the heat- 
transfer rate near the stagnation point approaches zero for an m Q of 0.4. Downstream, 
however, the effect of injection rate on heat-transfer rate decreases rapidly. Figure 37 
shows that the same conclusions apply to Stanton number distributions. 

Figure 38 shows the effect of mass injection on skin- friction coefficient distributions. 
In the mass- injection region, large reductions in skin friction occur. Downstream, how- 
ever, the skin-friction values with injection are approaching the no- injection values. For 
example, for a stagnation mass-injection rate of 0.4, the skin friction at % = 0.4 is 
90 percent less than the no-injection value, but at £ = 3.0 the reduction is only 9 percent. 

The effect of mass injection on wall pressure distribution is shown in figure 39. An 
almost negligible effect is obtained. The wall pressure is less with mass injection but 
the pressure differences never exceed 4 percent for the mass- injection rates considered. 
This result is due primarily to the way that the injected mass alters the shock shape and, 
consequently, the shock pressure. This condition is evident in figure 35, where the shock 
slope decreases with increasing injection rate. Hence, the local shock angle a (fig. 1) is 
less and so is the shock pressure. Obviously, the way mass injection alters the shock 
shape and shock pressure will depend on the assumed mass injection distribution. 

For some hypersonic flow conditions, in particular, hypervelocity planetary entries 
where radiant heating is very large, the mass-injection rates can be large enough to 
insure that the convective heating is reduced to zero. Calculations were made to demon- 
strate the present ability to calculate the effect of mass injection on stagnation heat trans- 
fer and shock shape at such injection rates. Figures 40 to 42 show such an effect for the 
injection of equilibrium air into reacting equilibrium air. Figure 40 shows the manner in 
which the stagnation wall heat-transfer rate decreases with increasing mass injection and 
becomes zero at a mass-injection rate of about 0.4. Figure 41 shows the influence of 
injection on stagnation shock standoff distance. For the range of injection rates shown, 
the shock standoff distance increases linearly with increasing injection rate. Figure 42 
shows the effect of even larger injection rates on shock standoff distance. (The results 
for injection rates greater than 0.4 did not consider the downstream influence on the stag- 
nation solution.) For injection rates greater than 0.4, the shock standoff distance continues 
to increase with increasing injection rate but the slope of the curve decreases. 
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CONCLUSIONS 


Equations describing reacting and nonreacting viscous shock layers about an axisym- 
metric body at an angle of attack of 0° are presented. These equations account for mass 
injection and multicomponent diffusion. An implicit finite-difference technique for solving 
these equations is discussed. Comparisons of the present results with previously reported 
solutions are made for a range of flow conditions that include shock Reynolds numbers 
(p*h* a * /Mg ) of 100 to 17 931. Calculated results are described which show the effect of 
chemistry model, surface catalysis, diffusion model, and mass injection on flow param- 
eters and surface transport. Some of the solutions provide information previously unavail- 
able for reacting air. For example, the downstream viscous -shock-layer solutions that 
treat the diffusion as multicomponent, with or without injection of air, were previously 
unavailable. These results were obtained for both nonequilibrium and equilibrium flow 
over slender bodies. 

Results of the study lead to the following conclusions: 

1. The chemistry model substantially influences flow parameters and surface trans- 
port. Results show that the frozen and equilibrium chemistry models overpredict wall 
heat-transfer rates for the flow conditions considered in this study. For example, the 
stagnation heating rate obtained with an equilibrium calculation is 54 percent greater than 
comparable results with a nonequilibrium and noncatalytic wall calculation. 

2. Wall catalysis significantly influences wall heat- transfer rate and species 
composition near the wall but has little influence on the other flow parameters considered. 
For an equilibrium catalytic wall, the wall heat-transfer rate for nonequilibrium flow is 
shown to be about the same as that for equilibrium flow. 

3. For reacting air, results for multicomponent and binary diffusion models were 
compared. The binary approximations included both constant Lewis numbers (1.0 and 
1.4) and a variable Lewis number, where the diffusion coefficient was the diffusion coeffi- 
cient for atomic oxygen diffusing into molecular nitrogen. Results show negligible differ- 
ences in most flow parameters and surface transport when binary diffusion is used in lieu 
of multicomponent diffusion. The largest effect of diffusion model is on heat-transfer 
rate and species distributions. When constant Lewis numbers of 1.0 and 1.4 are used, 
the heat transfer is within 10 percent of the results with multicomponent diffusion. When 
a variable Lewis number is used, the differences are negligible. This result applies to 
both nonequilibrium and equilibrium air. For multicomponent diffusion a small amount of 
elemental diffusional separation is observed, and consequently, a small difference in 
species profiles occurs. 

4. Mass injection significantly alters all flow parameters and surface transport. 

The effect that mass injection has on the flow field downstream of the injection region 
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decreases very rapidly with downstream distance. The solutions with mass injection 
clearly demonstrate the capability of the present analysis to study problems with 
distributed massive blowing. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 7, 1973. 
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APPENDIX A 


FACTORS USED TO NONDIMENSIONALIZE THE 
VISCOUS- SHOCK- LAYER EQUATIONS 

The viscous-shock-layer equations are nondimensionalized by using the following 
relations : 


* * 
u = uU 


* *2 
h = hU>_r 


* * 
v = vU 


. * * 
* W.PooUoo 
W- = — — g 

1 a 


* TU„ 
T =~^- 


r * i M ref 


* tt *2 


P* = P PooU, 


s* = sa 


* * 
p = PP m 


* * 
n = na 


p = mm„ 


* * 
K - KB 


* * * 
K = Kju -C 

ref p,° 


r* = ra 


where 


* * 

c = c c 

p p p,° 


Also, four dimensionless parameters appear in the shock-layer equations. They are 


* * 
C rj P 

TM = P 
^Pr r- 
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APPENDIX B 


THERMODYNAMIC AND TRANSPORT PROPERTIES 

Thermodynamic properties, Cp ? i and h*, and transport properties, /j* 5 K*, and J9*j, 
are required for each species considered. For equilibrium calculations, the free 
energies F i are also required. Since the multicomponent gas mixtures are considered to 
be mixtures of thermally perfect gases, the thermodynamic and transport properties for 
each species are calculated by using the local static temperature. 

Then the properties for the gas mixtures are determined in terms of the individual 
species properties. 


Thermodynamic Properties 

Values for the thermodynamic properties as a function of temperature are obtained 
by using polynomial curve fits for each chemical species. The following polynomial 
equations are used: 

Specific heat: 


= a x + a 2 T* + a 3 T* 2 + a 4 T* 3 + a 5 T* 4 
R* 


Enthalpy: 


* * *9 *9 *4 

h« agT a^T a^T a^T ^0 

- - _ = a* + — + — — + + — + _ 

R*T* 2 3 4 5 T* 


and 


(Bl) 


(B2) 


Free energy: 


*o _ 

F- * a«T 

= a, (1 - log e T ).h t 

R T 2 


*2 *2 _*4 
a 3 T . a 4 T _ a 5 T + a 6 " a 7 
6 12 20 T * 


(B3) 


where F. is the free energy of species i at 1 atmosphere pressure (standard state). 

The development of these curve fits and a tabulation of the polynomial constants (a^ to a^) 
are presented in reference 32. For the species used in this study, the constants are 
given in table IX. 
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APPENDIX B - Continued 
Transport Properties 

The gas mixtures considered are assumed to be mixtures of thermally perfect gases. 
The expressions used for calculating the transport properties of the pure species will be 
presented first, and are followed by the expressions used for calculating the transport 
properties of the mixture. 

To a first approximation (first-order kinetic theory), the viscosity for a perfect gas 
(ref. 33) is given as 



2.6693 x 10 


yf r 



a *2 n *( 2 , 2 ) 

i ij 


and the thermal conductivity for a monatomic gas is given as 


(B4) 


or 


* 

K 

i,mon 


1.9891 x 10" 4 yr^ /M* 


G. Cl 
l 


*( 2 , 2 ) 

ij 


* 

i/mon 



* 


R 


(B5) 


(B6) 


The thermal conductivity for a polyatomic gas has an additional contribution due to the 
transfer of energy between translational and internal degrees of freedom. Consequently, 
the expression for the thermal conductivity of polyatomic molecules is given (ref. 34) as 




M + 1.32 
4 




In the foregoing equations 
T* temperature, K 

M? molecular weight, g/mole 


(B7) 
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* 
0 , 

1 

collision diameter, A 

* 

M i 

viscosity, g/cm-sec 

* 

K i 

thermal conductivity, cal/ sec-K 

n*.< 2 ’ 2 > 

collision integral for viscosity 


The derivations of equations (B4) and (B5) are based on a very simple gas model (ref. 34); 
that is, a nonreacting mixture of monatomic gases at relatively low temperatures. How- 
ever, equations (B4) and (B7) have been found to be reasonably accurate for temperatures 
that are less than that required for ionization. Figures 43 and 44 show the viscosity and 
thermal conductivity values, respectively, for each chemical species. 

Expressions for the multicomponent viscosity and thermal conductivity according to 
rigorous kinetic theory are given in reference 33, However, these expressions are 
cumbersome and computationally time consuming. Consequently, most studies use 
approximate techniques for evaluating the mixture viscosity and thermal conductivity. 

The mixture viscosity is obtained by using the se.miempiri.cal formula of Wilke 
(from ref. 20) 


where 


N s 

i=l 


x i 


N 


E x i *ii 

3=1 



(B8) 


(B9) 


Figure 45 shows the mixture viscosity of equilibrium air as a function of temperature and 
for a pressure of 1.0 atmosphere (1 atmosphere = 101.3 kN/m^). Also shown is the 
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viscosity as predicted by Sutherland's equation. Note that the values of viscosity used in 
this study are greater than those predicted by Sutherland’s equation. For most of the 
calculations made in this study, the wall temperature is 1500 K. At this temperature, the 
viscosity predicted by the Sutherland equation is 7 percent less than the value used 
herein. 

The mixture thermal conductivity is obtained by a method analogous to that used for 
viscosity. The relation used (ref. 20) for calculating the mixture conductivity is 


K 


N * 

s x.K . 

* ii 


- L 


N 

i=l s 


r x. ■*.. 

t—i i ii 

j=i 


(BIO) 


where the 0.. coefficients are identical with those that appear in the viscosity equation. 

Values for the equilibrium air mixture thermal conductivity as a function of tem- 
perature and for a pressure of 1 atmosphere are given in figure 46. These values are also 
compared with the mixture conductivity obtained by assuming a constant Prandtl number 
equal to 0.72, Sutherland' s value of viscosity, and the specific heat as given by equa- 
tions (13) and (Bl). The resulting values of conductivity are less (10 percent less at 
1500 K) than the conductivity values used herein. 

The binary diffusion coefficients are obtained by using the relation given in 
reference 33 


(2.628 x 10 -3 )T* 3/2 


6 * = 

il 


* * 

M. + M. 

“ * *~ 

2M t M. 


(Bll) 


pV^M) 

1] 1J 


where .£2^ 

a *(2,2) 11 




is the collision integral for diffusion. The collision integrals £2, 


‘( 1 , 1 ) 


ij 


are functions of the nondimensional reduced temperature 


il 


and 


T + 



(B12) 
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* j(c 

where e.. is the maximum energy of attraction between colliding molecules and k is 

* * 

Boltzmann’s constant. The parameters and e.j are estimated satisfactorily by 


a* = I (a* + a*) 

ij 2 1 3 


and 


* / W'"* 

e. . = / e. e. 
V i 3 


(B13) 


(B14) 


The force constants o j and e^/k for each species considered in this study are given 
in table X. The force constants are those given in reference 34. The collision integrals 
were obtained from reference 33 (pp. 1126-1127) where they are tabulated as a function of 
reduced temperature T + . These collision integrals are based on the Lennard Jones 6-12 
potential for the interaction of colliding molecules. 

For the binary diffusion calculations, the diffusion is approximated by two methods. 
One method assumes that all species have diffusion coefficients equal to an effective binary 
coefficient. The effective binary coefficients are determined for a preassigned diffusing 
pair of species. For example, the effective binary diffusion coefficient for air will be that 
of atomic oxygen diffusing into molecular nitrogen evaluated at the local temperature and 
pressure according to equation (Bll). The second binary diffusion approximation is that 
the Lewis number is some constant value. Both approximations are used in this study. 

The multicomponent diffusion coefficients are dependent on the concentration of the 
species. These coefficients are obtained by using the relations given in reference 1, which 
may be written as 




% 


* 

M. _* 

-Ki 


(B15) 


where the quantities are coefficients in the matrix which is the inverse of the matrix 
with the following coefficients: 


C. 


i+ M* 
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K « 
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jMI ' Me Sie 


(i ± j) 


(i = j) 


(B16) 
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APPENDIX C 


TRANSFORMED VISCOUS-SHOCK- LAYER EQUATIONS 

This appendix presents the transformed viscous-shock-layer equations and boundary 
conditions. First, the relations defining the transformed variables and coordinates are 
given. Next, the general equations and boundary conditions are given. Then the special 
form of the equations for the stagnation streamline are developed along with the stagnation 
shock relations. 


Transformation Relations 


To simplify the numerical computations, a transformation is applied to the viscous- 
shock-layer equations. This transformation is accomplished by normalizing most of the 
variables with their local shock values. When the normal coordinate is normalized with 
respect to the local shock standoff distance, a constant number of finite-difference grid 
points between the body and shock are used. This procedure eliminates the need for inter- 
polating to determine shock shape and the addition of grid points in the normal direction 
as the computation moves downstream. 

The transformed variables are 





The transformations relating the differential quantities are 
, n* , 

3 _ 3 S 3 

3s 3£ n g dr] 


(C2) 
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where 


n s = 


dn 


d{ 


_a_ _ 1 

an n 


s 


a 


(C3) 

(C4) 


and 


a 2 = 1 


an n„ 


a v 


(C5) 


11 11 

The transformations used to express the Shock- oriented velocities u g and v g in terms 
of the body- oriented coordinate system (fig. 1) are 


u = u sin (a + p) + v cos (a + p) 

SS S 


(C6) 


and 


v o = -u cos (a + p) + v sin (a + P) 
s s s 


(C7) 


Transformed Equations 

After the governing equations are written in transformed variables and coordinates, 
the second- order partial differential equations are written in the following form: 


a 2 W 

an 2 


aW 


aW 


+ fll ^W + a 3 + a 4 - =0 


3f 


(C8) 


The quantity W represents u in the s-momentum equation, T in the temperature energy 
equation, H in the enthalpy energy equation, in the species continuity equations, and C g 
in the elemental continuity equation. The coefficients to are written as follows: 
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s- momentum, W = u: 
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The preceding energy equation is for the thin shock-layer approximation. When equa- 
tion (3) is used for the n- momentum equation, the following term must be added to 
(eq. (Cllc)) 
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Species continuity, W = C.: 
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where 


and 


v= - 


jM 

e 2 PL i 


(C12b) 


a 3 " 


PL i 


3 PM. 


1 + n PM. 


j cos d 


. o 2 
w i n s 


B T] 


S 1 


1 + n t? k r + n tj cos 8 2 

\ S S / v 


(C12c) 


a 4 = - 


n s p s U s PU 


e^PL t (1 + n V k) 


(C12d) 


PL. = 


^ N Le,i 

N Pr,s *Pr 


(C12e) 


N 


PM t 


M g M 


3 C q 


— __ £ Ab ik — - 
N Pr,s N Pr k=l dr) 

Mi 


(C12f) 


for multicomponent diffusion and 


PL. = PL = 
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and 
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for binary diffusion. 
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Elemental continuity, W = : 
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for binary diffusion. 
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The remaining equations are written as follows: 
Global continuity: 


— jn g (r + n g n cos o9 P g u g pu] + JL |(r + n g n cos 6$ [(1 + n g w) p g v g pv 

H V p s u s pl ]} = 0 (C 14 ) 

n- momentum: 
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if the thin- shock layer approximation is made, and 
State: 
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(C15a) 


(C15b) 


(C16) 


The terms Wj and Wg which appear in the temperature form of the energy equation 
and the terms w? and w| which appear in the species continuity equations are quantities 
that involve the rate of production terms, w.. As discussed in references 1, 7, and 11, the 
way the production terms are written is very important in achieving convergence of the 
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iteration procedure. Consequently, for the energy equation, the production terms are 
written as (see ref. 7) 




'w.' 


k+1 


p/ k \9T P/ k 


t). ( T k + 1 ■ T k) 


(C17) 


where k denotes the iteration for which the solution is known and k+1 the iteration for 

N 

s 

which a solution is required. The term ^ h.w. which appears in equation (4) is written 

1=1 

as follows: 
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E Vl = *81*1 + T s w s2 


(C18) 


i=l 

and appears in equations (ClOb) and (ClOc). As for the species continuity equations, the 
production term is written as 


w i 


= w. 


° - C A 


(C19) 


and appears in equations (C12b) and (C12c). Hence, equations (C17) and (C19) express the 
production terms in terms of temperature for the energy equation and in terms of species 
mass fraction for the species equations. Blottner in reference 1 comments on the neces- 
sity of using expressions like equations (C17) and (C19). 

For frozen, equilibrium, and nonequilibrium chemistry, equations (C8) to (C16), 
along with the appropriate boundary conditions and relations for the thermodynamic and 
transport properties, are the governing relations used to describe the viscous shock layer. 


Boundary Conditions 

Conditions at the body surface .- The surface boundary conditions in terms of trans- 
formed variables are 
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The surface species boundary conditions (eq. (C25)) are written as 
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for multicomponent diffusion and 
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for binary diffusion. 

The elemental boundary conditions (eq. (B26)) are written as 
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f'J r\, 


(C36) 


for binary diffusion. When the diffusion is binary and the injectant has the same composi- 
tion as the free stream, the elemental concentrations are constant and equal to the free- 
stream value. For this case, the elemental equations are superfluous and the surface 
conditions become 



(C37) 


However, when the diffusion is multicomponent, the elemental distribution along the normal 
to the surface is not constant and the elemental conservation equations must be used to cal- 
culate the elemental distribution. 

Conditions at the shock. - The shock conditions are determined by solving equa- 
tions (35) to (42). The transformed shock conditions become 


u = T = H = v = p = p = 1 


(C38) 


at i? = l. Note that the elemental species mass fractions are not normalized with their 
respective shock values. 


Stagnation Streamline Equations 

When downstream numerical solutions are required, it is necessary to have an 
accurate solution for the flow along the stagnation streamline. A truncated series, which 
has the same form as that used by Kao in reference 35, is used to develop the stagnation 
streamline equations. The flow variables are expanded about the axis of symmetry with 
respect to the nondimensional distance % along the body as follows : 


p(f ,n) - p 1 (n) + p 2 (n ) r 2 + • • • 

(C39a) 

u(£,t?) = u 1 (n)? + • • • 

(C39b) 

v(£ ,i?) = Vj(ti) + . , . 

(C39c) 

P U , v) = p 1 (t?) + • • • 

(C39d) 
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T(f,i?) = T 1 (tj)+ . . . (C39e) 

h( I, tj) = hj( ^ + • • • (C39f) 

M (I ) = M^(i? ) + • • • (C39g) 

K(|,tj)= • • • (C39h) 

C p (M) = C pjl (i?)+ ’ * • (C39i) 

G s (l,r,)=^ 1 j(n)+..- (C39j) 

C i U,r ? )=C i>1 (T 1 ) + ... (C39k) 

The shock standoff distance is written as 

n s = n ls + n 2s £ 2 + * * * (C40) 


Furthermore, £ is small and the curvature k is approximately one in the stagnation 
region. Consequently, the geometric relations (see fig. 1), including terms of order £, can 
be written as 

P « £ (C41) 

and 


a 



2n, 


2 s 


1 + n- 


Is 


(C42) 


Therefore, 


sin (a + p ) * 1 


(C43) 
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APPENDIX C - Continued 


and 


2n 2 J 

COS (a + |3 )<=« _ — 


1 + n u 


(C44) 


The shock relations (eqs. (35) to (42)) in terms of expanded variables become 


v s = v ls + " * 


Is 


(C45) 


and 


u = u, % + 
s Is 


2n, 


2s ( 1 _± 


1 + n ls \ P ls/ 


(C46) 


p s = p ls + p 2s f2 + 11 

7 Mi 


Is/ 


- r 


i - 


2n, 


2 s 


Is/ 


1 + n 


Is/ 


(C47) 


h s - h ls 


+ I n 


M? (* - 1) 2 


Is, 


(C48) 


Equations (40) to (42) are unchanged. 

By examining these equations, it is seen that the equations for u g and p g contain 
ngg. This term cannot be determined from the stagnation solutions since it is a function 
of the downstream flow. Consequently, a value must be assumed for ngg. In this study, it 
was assumed to be zero to start the solution, but this assumption is then removed by iter- 
ating on the solution by using the previous shock standoff distances to define n 2s . The 
effect of the downstream shock shape on the stagnation-point solution is elliptic rather 
than parabolic. 

Along the stagnation streamline, the second- order differential equations are written 
as 
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APPENDIX C - Continued 



+ 


u w 

a * + a 0 w + a„ = 0 

1 d 2 3 


The coefficients are defined as 
s-momentum; W = u: 


- 1 d "l “is “Is "ls T !s "l v l 


a l - 


^ di? 


+ (j + 1) 


1 + n lB , 


e /z 


Is 


n 


a 2 “ 


Is 


1 + 


n ls* 


i 1*L + 

^ d») 1 + n^ g n 


(j + 1) + 


p ls n ls u ls P 1 P 1 


e 2 M 


Is 


n ls p ls v ls p l v l 


• 2 "is U 1 


“3 ■ 


- 2p ls n ls 


e "is (1+ n ls’ l)u ls" 


- p 2 s p l n ls” dp l 

P 2 + - 

Pis n ls d7? 


Energy (temperature); W = T: 


a l = 


_1_ 

Kj drj 


+ (j + 1) 


n 


Is 




1 + n ls " K ls 


k L Wi 

1 i=l 


_ n ls P ls c p,ls v s P l c p,l v l 


e K 


Is 


K, 


(C49) 


(C50a) 


(C50b) 


(C50c) 


(C51a) 
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APPENDIX C - Continued 


a 2 = - 


n ls W s2 W 2 

e2K ls E l 


n l S v l S Pt s ? l d Pl 


n ls w sl 


eV K, K, dTJ 
Is Is 1 


e 2 T 


ls K ls K l 


w, 


Energy (enthalpy); W = H: 


(C51b) 


(C51c) 


= 1 d gl t dN pr>1 ^ (Hl)n ls n ls /> ls^Pr,ls v ls ? lV,l ? l 
dr? Np r j dr? 1+n^r? p 2 M-,., 


(C52a) 


Un^rJ 


« 2 M ls 


Mi 


a 2 = 0 


(C52b) 


_ N Pr,ls n ls N Pr,l 


M, H i 
Is Is 


Mi 


1 d^ + (j + 


n ls d,t 1+ ®1, 


(C52c) 


Species continuity; W = C.: 


dPL. 


ai =. 


PL. dr? 


1 + (j + 1) 


n 


Is 


1 + r?n 


Is 


^ls^lsVl 

e 2 PL. 


(C53a) 


a 2 


2 *1 
n ls w i 

e 2 PL. 


(C53b) 


a 3 


PL i 


dPM. 

dr? 


- + (]'+ 1) 


PM i n ls 


.o 2 " 
w i n ls 


1 + r?n 


Is 


(C53c) 
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APPENDIX C - Continued 


Elemental continuity; W = C s : 


a l 


J_ dPL 
PL 


+ (J + D 


Is 

1 + tj n 


Is 


n ls p ls v ls p l v l 

9 (M 

t PL 


a 2 “ 0 


CLq “ 


PL 


dPM 0 n 1o PM„ 

' + a + i) ls ‘ 


drj 


1 + n n 


Is 


The remaining equations are written as follows: 


Global continuity : 


A [» + n ls „ )i +1 o lg v lg ^vJ 

= -(j + l)n ls (l + n ls T,)! 

n- momentum: 

d Pl v ls p ls - - *1 

= - = PiV, 

dT? P ls drj 


(C54a) 

(C54b) 


(C54c) 


(C55) 


(C56a) 


When the thin shock- layer approximation is made, the n- momentum equation becomes 


_- 0 

d 77 

The Pg term that appears in equation (C50c) can be expressed as 


(C56b) 
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APPENDIX C - Concluded 


dp 2 = 

p ls u ls n ls 

p l u l 

^ 2p ls u ls n 2s v ls 


dVj 

dr? 

Pis 

1 + 1 “is 

Pis 

» ' "“is 

dv 



dp 2 

drj 


A ls u l,s n ls 



1 + 



(C56d) 


These equations along with the equation of state constitute the nonlinear ordinary differen- 
tial equations that are solved along the stagnation streamline. 
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TABLE I.- THIRD BODY EFFICIENCIES RELATIVE TO ARGON 


Catalytic 

bodies 

Z (J-N s ),i 

Mi 

M 

m 2 

2,i 

m 3 

3,i 


Efficiencies relative to argon of — 


0 

°2 

N 

N 2 

NO 

Q 

i = 2 

i = 3 

i = 4 

1 = 5 

25 

9 

1 

2 

1 

1 

1 

0 

2.5 

1 

20 

1 

20 

20 

20 


TABLE II.- REACTION RATE COEFFICIENTS 


Reaction, 

r 


l/’mole\" Q!r 

s \cm 3 / 

3.61 x 10 18 
1.92 x 10 17 
4.15 x 10 22 
3.97 X 10 20 
3.18 x 10 9 
6.75 X 10 13 


59.4 
113.1 
113.1 

75.6 

19.7 

37.5 


1 An ol e V 
s \cm 2 J 


3.01 x 10 15 
1.09 X 10!6 
2.32 X 10 21 

1.01 X 10 20 
9.63 x 10 11 
1.50 X 10 13 































TABLE IB.- ANALYSES FOR WHICH COMPARISONS ARE MADE 


Investigator(s) 

Analysis 

Air chemistry 
model 

Diffusion 

Region of 
application 

Whitehead, 
reference 8 

Viscous shock 
layer 

Frozen and 
equilibrium 

Binary 

Stagnation and 
downstream 

Edelman and Hoffman, 
reference 15 

Viscous shock 
layer 

Equilibrium 

Binary 

Stagnation 

Goldberg and Scala, 
reference 16 

Viscous shock 
layer 

Equilibrium 

Binary 

Stagnation 

Anderson and Lewis, 
reference 28 

Boundary layer 

Equilibrium 

Binary 

Stagnation 

Blottner, as given 
in reference 28 

Boundary layer 

Nonequilibrium 

Multicomponent 

Stagnation and 
downstream 

Smith, as given 
in reference 28 

Boundary layer 

Nonequilibrium 

Binary 

Stagnation and 
downstream 

Blottner, 
reference 11 

Viscous shock 
layer 

Nonequilibrium 

Binary 

Stagnation 


TABLE IV.- FREE-STREAM CONDITIONS 


Altitude, 

km 

Velocity, 

km/s 

Temperature, 

K 

Pressure, 

atm 

Density, 

g/m 3 

Mach 

number 

Free -stream 
Reynolds 
number 
(a) 

85.34 

7.92 

180.65 

3.82 x 10" 6 

7.46 X lO' 6 

29.40 

1.482 x 10 3 

76.20 

6.10 

195.4 

2.01 x 10" 5 

3.62 x 10" 5 

21.75 

5.170 X 10 3 

60.96 

6.10 

253.88 

1.95 x 10 “ 4 

2.70 X 10“ 4 

19.08 

3.116 X 10 4 

45.72 

6.86 

266.67 

1.28 x 10 "3 

1.70 x 10‘ 3 

20.98 

2.217 X 10 5 

30.48 

6.10 

266.98 

1.09 x 10 -2 

1.70 X lO" 2 

20.18 

2.149 X 10 6 


a Based on a 0.3048 -meter nose radius. 
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TABLE V.- FROZEN AIR STAGNATION SHOCK AND WALL VALUES 


Altitude, 60.960 km; 
U * = 6.096 km/s; 
a* = 2.54 cm; 

NRe, co = 2.597 X 10 s ; 
_N Re ,s = 1-914 x 102_ 


Properties 

Value at shock 

Value at wall for - 

n s = 9.079 x 10" 2 ; 
q* = 4.865 MW/m 2 

n s = 9.215 X 10 -2 ; 
q* = 4.796 MW/m 2 

Temperature, K 

1.331 x 10 4 

1.000 X 10 3 

1.500 X 10 3 

Pressure, atm 

8.822 X IQ" 2 

9.439 x 10-2 

9.439 X 10-2 

Density, kg/m 3 ...... , . . 

2.331 x 10-3 

3.317 x 10-2 

2.211 X 10-2 

Viscosity, N-s/m 2 . 

2.197 x 10" 4 

4.077 X 10' 5 

5.294 X 10-5 

Velocity, km/s 

-7.096 x 10“1 

0.0 

0.0 

Enthalpy, MJ/kg 

1.857 X 10l 

7.630 X 10“1 

1.355 x 10° 

Molecular weight, g/g-mole . . 

2.888 x 10l 

2,888 x 10l 

2.888 X 10l 

Prandtl number 

7.018 x 10-1 

6.938 x 10-1 

6.961 X 10-1 

Mass fractions: 




02 

2.400 x 10-1 

2.400 x 10-1 

2.400 X 10"! 

n 2 

7.600 x 10-1 

7.600 x 10"1 

7.600 X lO'l 
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TABLE VI.- EQUILIBRIUM AIR STAGNATION SHOCK 
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TABLE VI.- EQUILIBRIUM AIR STAGNATION SHOCK 
AND WALL VALUES - Continued 

(b) Altitude, 85.344 km; 

U* = 7.925 km/s; 
a* = 30.480 cm; 

N Re ,co = 1-482 XlO 3 ; 

N Re> s = 1-498 Xl02; 

n s = 3.801 x lO" 2 ; 

q* = 7.861 X 10 _1 MW/m 2 


Properties 

Value at 
shock 

Value at 
wall 

Temperature, K . 

•5.733 

X 

10 3 

1.349 

X 

10 3 

Pressure, atm V 

4.396 

X 

10 -3 

4.528 

X 

10" 3 

Density, kg/m 2 

1.485 

X 

10“ 4 

1.186 

X 

10' 3 

Viscosity, N-s/m 2 . 

1.203 

X 

10" 4 

4.940 

X 

10" 5 

Velocity, km/s . . . , 

43.981 

X 

10“ 1 

0.0 



Enthalpy, MJ/kg 

3.163 

X 

ioi 

1.176 

X 

ioo 

Molecular weight, g/g-mole 

1.580 

X 

10 1 

2.888 

X 

10 1 

Lewis number 

1.400 

X 

10° 

1.400 

X 

10° 

Prandtl number \ 

6.734 

X 

10“1 

6.957 

X 

10" 1 

Mass fractions: 

} 

\ 





10" 8 

0 . . . 

2.398 

X 

10 _1 

5.331 

X 

°2 • • • • 

2.098 

X 

10" 6 

2.396 

X 

10' 1 

N 

5.927 

X 

io-i 

5.293 

X 

lo - 13 

n 2 

J..671 

X 

10 _1 

7.597 

X 

10" 1 

NO 

^.776 

X 

10“ 4 

6.008 

X 

10‘ 4 
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TABLE VI.- EQUILIBRIUM AIR STAGNATION SHOCK 
AND WALL VALUES - Continued 


(c) Altitude, 45.720 km; 

U* = 6.864 km/s; 
a* = 1.27 cm; 

N Re ,oo = 9.236 x 10 3 ; 
Nrs, J- 1-134 X 103; 
n s = 5.601 x . 10 "2j 
q* = 2.695 x 10 1 MW/m2 


Properties 

Value at 
shock 

Value at 
wall 

Temperature, K ................. 

6.783 x 10 3 

1.476 x 10 3 

Pressure, atm . . .... . . 

7.729 x 10" 1 

8.038 x 10' 1 

Density, kg/m 3 ........ .I....... . 

2.592 x 10-2 

1.916 x 10' 1 

Viscosity, N-s/m^ ...... . | . 

1.369 X 10" 4 

5.237 x 10-5 

Velocity, km/s 

-4.715 x 10" 1 

0.0 

Enthalpy, MJ/kg 

2.370 X 10 1 

1.258 x 10° 

Molecular weight, g/g-mole ........... 

1.863 X 10 1 

2.888 x 10l 

Lewis number . . . . 

1.400 x 10° 

1.400 x 10° 

Prandtl number ................. 

6.627 x 10 _1 

6.958 x 10" 1 

Mass fractions: 



0 ....... 

2.371 x 10 _1 

3.638 x 10 -7 

02 * 

1.108 x 10“ 4 

2.397 x 10 _1 

N 

3.260 X KT 1 

3.983 x 10 -15 

n 2 . . 

4.310 x 10 _1 

7.593 X 10' 1 

NO ... . . . . . . . 

5.712 x 10' 3 

9.205 X 10“ 4 
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TABLE VI.- EQUILIBRIUM AIR STAGNATION SHOCK 
AND WALL VALUES - Continued 


(d) Altitude, 76.20 km; 

U* = 6.096 km/s; 
a* = 5.532 cm; 

NRe,« = 9-383 X 10 2 - 
N Re,s = '-O' 50 * !0 2 ; 
n s = 5.298 x 10“ 2 ; 
q* = 1.475 X 10° MW/m 2 


Properties 

Value at 
shock 

Value at 
wall 

■ ■ ■ - : 

Temperature, K 

5.291 x 10 3 

1.500 x 10 3 

Pressure, atm • 

1.245 x 10 -2 

1.290 X 10' 2 

Density, kg/m 2 

5.576 X 10" 4 

2.305 x 10' 3 

Viscosity, N-s/m 2 

1.221 X 10" 4 

6.434 X 10" 5 

Velocity, km/s 

-3.959 x 10 _1 

0.0 

Enthalpy, MJ/kg 

1.863 x 10 1 

1.938 X 10° 

Molecular weight, g/g- mole 

2.885 X 10 1 

1.970 X 10 1 

Lewis number 

1.000 x 10 1 

1.000 X 10 1 

Prandtl number 

7.019 x 10" 1 

6.815 X 10"! 

Mass fractions: 



0 

2.393 x 10" 1 

1.357 x 10 -3 

°2 

9.240 x 10" 6 

2.346 X 10 _1 

N 

2.423 X 10" 1 

1.490 X 10 -9 

n 2 

5.171 x 10-1 

7.565 X 10" 1 

NO 

1.267 X 10 -3 

7.535 X 10" 3 
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TABLE VI.- EQUILIBRIUM AIR STAGNATION SHOCK 
AND WALL VALUES - Concluded 


(e) Altitude, 30.48 km; 

U* = 6.096 km/s; 
a* = 2.54 cm; 

N R e,oo = 1-791 X10 5 ; 
N Re,s = i- 794 x io 4 ; 
n s = 6.755 X IO' 2 ; 
q* = 3.725 X 10 1 MW/m 2 


Properties 

Value at 
shock 

Value at 
wall 

Temperature, K 

6.967 X 10 3 

1.400 x 10 3 

Pressure, atm 

5.767 x 10° 

6.032 x 10° 

Density, kg/m 2 

2.084 x 10 _1 

1.516 x 10 1 

Viscosity, N-s/m 2 

1.475 x IO' 4 

5.191 x IO" 5 

Velocity, km/s 

-4.959 x IO' 1 

0.0 

Enthalpy, MJ/kg 

1.867 X 10* 

1.225 X 10° 

Molecular weight, g/g-mole 

2.897 X 10 1 

2.076 X 10 1 

Lewis number a 



Prandtl number 

6.835 X 10" 1 

6.992 X 10" 1 

Mass fractions: 



O 

2.325 x IO -1 

2.579 X IO" 8 

o 2 

3.571 x IO -4 

2.670 X 10" 1 

N 

1.761 X 10" 1 

1.105 X 10" 29 

n 2 

5.777 X 10' 1 

7.321 x 10" 1 

NO 

1.327 X 10" 2 

9.093 X 10" 4 


a Multicomponent diffusion. 
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TABLE VII.- NONEQUILIBRIUM AIR STAGNATION SHOCK AND WALL VALUES 
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a Lewis number was constant at 1.4. 




TABLE VHI.- NONEQUILIBRIUM AIR STAGNATION WALL CONDITIONS FOR A NONCAT ALYTIC WALL 
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a See table VII for shock properties.. 




TABLE DC.- CONSTANTS FOR POLYNOMIAL APPROXIMATIONS 
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(a) Mass balance at gas-solid interface. 
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(b) Energy balance at gas-solid interface. 
Figure 2.- Surface mass and energy balances. 
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(a) Shock in chemical equilibrium. 

Figure 4.- Flow chart for shock solution procedure. 
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(b) Shock frozen at free -stream chemistry. 
Figure 4.- Concluded. 




(a) Equilibrium chemistry. 

Figure 5.- Flow chart lor solution sequence of viscous -shock-layer equations. 
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(b) Frozen and nonequilibrium chemistry. 
Figure 5.- Concluded. 
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Present analysis for frozen air 
Whitehead (ref. 8) 







(c) Heat transfer variation with distance along body surface. 
Figure 6.- Concluded. 
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Present analysis for equilibrium air 
Whitehead (ref. 8) 


(a) Shock standoff variation with distance along body surface. 

Figure 7.- Comparison of present calculations with results of Whitehead for equilibrium 
air at an altitude of 60.96 km. U* = 6.10 km/s; = 1000 K; a* = 2.54 cm; binary 
diffusion; L = 1.0. 
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Present analysis 

Edelman and Hoffman (ref. 15) 



Temperature, K 

(a) Stagnation temperature profiles for an altitude 
of 85.34 km. u£ = 7.92 km/s; a* = 30.48 cm; 

= 1350 K; Lewis number = 1.4. 

Figure 8.- Comparison of present calculations with stagnation 
results of Edelman and Hoffman for equilibrium air. 



Mass fractions 



(b) Stagnation species profiles for an altitude 
of 85.34 km. U£> = 7.92 km/s; a* = 30.48 cm; 
T-^ = 1350 K; Lewis number = 1.4. 

Figure 8.- Continued. 
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Figure 9.- Comparison of present results with stagnation results of Goldberg and 

Scala for equilibrium air at an altitude of 76.20 km. u£> = 6.10 km/s; = 1944 K: 
Lewis number = 1.0; Nr 6 s = 100. 
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4 

Figure 10.- Skin-friction distributions for a 10° half-angle hyperboloid at an altitude 
of 30.48 km. U* = 6.096 km/s; = 1400 K; a* = 2.54 cm. 
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Temperature, K 

(a) Stagnation temperature profiles. 

Figure 11.- Comparison of present calculations with stagnation results of 

Blottner for nonequilibrium air at an altitude of 60.96 km. U^, = 6.10 km/s; 
T™ = 1000 K; a* = 2.54 cm; Lewis number = 1.4: Catalytic wall. 
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— — — Present analysis . 
Blottner (ref. 11) 
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(c) Stagnation species profiles. 
Figure 11.- Concluded. 
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Mass tractions, 



Figure 14.- Species profiles for nonequilibrium air with multicomponent diffusion 
and no mass injection at an altitude of 60.96 km. u£> = 6.10 km/s; = 1500 K; 
a* = 2.54 cm; noncatalytic wall. 
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(a) |=0.0. 


Figure 15.- Species profiles for equilibrium air with multicomponent diffusion and 
no mass injection at an altitude of 60.96 km. U£, = 6.10 km/s; = 1500 K; 
a* = 2.54 cm. 
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Figure 16.- Stagnation equilibrium and frozen air species and temperature profiles at 
an altitude of 60.96 km. U* = 6.10 km/s; T* = 1500 K; a* = 2.54 cm. 
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(a) | = 0.0. 

Figure 17.- Effect of chemistry model on tangential velocity profiles with 
multicomponent diffusion and no mass injection at an altitude of 60.96 km. 
U^j = 6.10 km/s; = 1500 K; a* = 2.54 cm; noncatalytic wall. 
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(a) | = 0.0. 4, 

Figure 18.- Effect of chemistry model on temperature profiles with multicomponent 
diffusion and no mass injection at an altitude of 60.96 km. u£, = 6.10 km/s; 

= 1500 K; a* = 2.54 cm: noncatalytic wall. 
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(b) £ = 2.0. 
Figure 18.- Concluded. 
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(a) | = 0.0. 

Figure 24.- Effect of wall catalyticity on species profiles at an altitude of 60.96 km. 
U* = 6.10 km/s; = 1500 K; a* = 2.54 cm; multicomponent diffusion. 
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Mass fractions 



(b ) | = 1.0. 
Figure 24.- Continued. 
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Nonequilibrium, noncatalytic wall 
Nonequilibrium, equilibrium catalytic wall 
Equilibrium, noncatalytic wall 


Figure 26.- Effect of wall catalyticity on nondimensional heat transfer at an altitude 
of 60.96 km. U* = 6.10 km/s; = 1500 IQ a* = 2.54 cm; multicomponent 
diffusion. 
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(a) £ = 0.0. 

Figure 29.- Effect of diffusion model on species profiles for equilibrium air 
at an altitude of 60.96 km. U£, = 6.10 km/s; = 1500 K; a* = 2.54 cm. 
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(a) |=0.0. 


Figure 30.- Effect of diffusion model on species profiles for nonequilibrium air 
at an altitude of 60.96 km. Uj> = 6.10 km/s; = 1500 K; a* = 2.54 cm. 
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(a) | = 0.0. 


Figure 32.- Comparison of nonequilibrium species profiles with and without mass 
injection of equilibrium air at an altitude of 60.96 km. U* = 6.10 km/s; 

= 1500 K; a* = 2.54 cm; multicomponent diffusion; noncatalytic wall. 








Figure 33.- Stagnation tangential velocity ratio profiles for different 
injection rates of equilibrium air into reacting nonequilibrium air 
at an altitude of 60.96 km. U* = 6.10 km/s; = 1500 K; 
a* = 2.54 cm; multicomponent diffusion; noncatalytic wall. 
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Figure 42.- Stagnation shock standoff distance as influenced by large mass injection 
rates of equilibrium air into reacting equilibrium air at an altitude of 60.96 km. 
U£> = 6.10 km/s; = 1500 K; a* = 2.54 cm; binary diffusion. 
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Figure 45.- Viscosity for equilibrium air at a pressure of 1 atmosphere. 




Figure 46.- Thermal conductivity of equilibrium air at a pressure of 1 atmosphere. 
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